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Abstract 

^^ . The stationary one dimensional Schrodinger-Poisson system on a bounded in- 

^>| terval is considered in the hmit of a small Debye length (or small temperature). 

Electrons are supposed to be in a mixed state with the Boltzmann statistics. Using 

various reformulations of the system as convex minimization problems, we show that 
C^^ , only the first energy level is asymptotically occupied. The electrostatic potential is 

shown to converge towards a boundary layer potential with a profile computed by 

means of a half space Schrodinger-Poisson system. 
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> : 1 Introduction and main results 

00 

2^ ■ 1.1 Introduction 

VQ ■ The Schrodinger-Poisson system is one of the most used models for quantum transport 

of charged particles in semiconductors as well as for quantum chemistry problems [Sj El 
[Hl[ISllIl[I51[ISllIHll2ni[2Ill22l[2SlEn]. it describes the quantum motion of an ensemble 
of electrons submitted to and interacting with an electrostatic potential. The electron 
ensemble might be completely confined or in interaction with reservoirs. In the latter 
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Vh ' case, one speaks about open systems for which the particles are described by means of 

the scattering states of the Schrodinger Hamiltonian corresponding to the electrostatic 
potential which is in turn coupled to electron particle density through the electrostatic 
interaction. This leads to nonlinear partial differential equations whose analysis involves 
scattering theory techniques and limiting absorption theorems [231 IH E] and in which the 
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repulsive character of the electrostatic interaction plays an important role in the analysis 
(it provides the necessary a priori estimates for solving the problem). 

For closed systems, the particles are described thanks to the eigenstates and eigenen- 
ergies of the Schrodinger Hamiltonian. The electron density is the superposition of the 
densities of the eigenstates with an occupation number decreasingly depending on their 
eigenenergy. The coupling is again obtained through the Poisson equation modeling the 
electrostatic interaction. This problem was reformulated by Nier [20| [2T| [22] as a mini- 
mization of a convex function (whose unknown is the electrostatic potential) which allows 
us to prove existence and uniqueness results. In [15j, one can find generalizations in- 
cluding local contributions to the potential and which can be included in the functional 
to be minimized. This short review partially covers stationary problems. For evolution 
problems, an extended bibliography is available, and we refer the reader to the books of 
Markowich, Ringhofer, and Schmeiser [TH] and Cazenave PH] for references. 

In this paper we are interested in a singularly perturbed version of the Schrodinger- 
Poisson system which arises from the description of the so-called two dimensional electron 
gases [UllI]. The electrons, in such systems, are strongly confined in one direction, at the 
interface between two material, and are free to move in the two remaining ones. In [6], 
the analysis of the Schrodinger equation of strongly confined electrons in one direction 
is performed. The confined direction is called z and the confining potential is assumed 
to be given and scaled as ^V^(-), where e is a small parameter. Approximate models 
for the transport direction (orthogonal to z) derived heuristically in the previous works 
[25| |26| |27] are then analyzed in [6l |2l]. The aim of the present work is to somehow justify 
the scaling 4-K(-) by the analysis of the self-consistent Schrodinger-Poisson system in 
the z direction. This is why we shall forget about the transport issues in the orthogonal 
direction and assume that the considered system is invariant with respect to it. The 
parameter e in the present work is linked to the scaled Debye length as shall be explained 
later. The analysis relies on the minimization formulation of the problem leading to a 
singularly perturbed functional. After a rescaling argument, we are led to the analysis 
of a half space Schrodinger-Poisson system in which only the first eigenstate is occupied. 
Additional estimates are obtained thanks to reformulation of the single state Schrodinger- 
Poisson system as another minimization problem whose unknown is the first eigenfunction 
(and not the potential). This formulation is used in quantum chemistry [9]. 

Let us now come to the precise description of the problem and the results. The system 
is one dimensional and occupies the interval [0,1]. The electrostatic energy is given by 
14(2;). It satisfies the following one dimensional stationary Schrodinger-Poisson system: 

— -TY ^y^v = ^v'^v^ 2;e[0,l], 

V?p G E^{^, 1), (/?p(0) = 0, (^p(l) =0, / Lp^Lpq = Spg, 

p=i p=i 

y(o) = o, ^(i) = o. 

dz 

The dimensionless parameter e is a small parameter which is devoted to tending to zero. 
The choice of the third power is done for notational convenience as shall be understood 



later. This parameter is related to the Debye length and shall be explicitly given by the 
rescaling of the Schrodinger-Poisson system ( 1191) (see subsection 1.3). The eigenvalues 
of the Schrodinger operator {Sp)p are the energy levels in the potential well. The sum 
in the right-hand side of the Poisson equation includes all eigenvalues of the Schrodinger 
operator. In the limit e — )■ 0, one expects that the wave functions concentrate at 2 = 0. 
The boundary condition for the potential at 2: = 1 is physically justified in some physical 
situations such as in bulk materials. However, a Dirichlet condition is more commonly 
used in such problems. The analysis can be carried out in that case with the cost of 
technical complexity since a new boundary layer at 2; = 1 will appear and the eigenvalues 
will have asymptotically a double multiplicity. For simplicity, we do not consider this 
case. Since the density is very high in the limit e — > O"*", the Boltzmann statistics should 
be replaced by the Fermi-Dirac ones. The analysis can be done in this case with the 
cost of technical complications. More detailed comments about this are given in the last 
section of this paper. In order to analyze the boundary layer, we make the change of 
variables 
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Then, U verifies — ^^ = -j Ylt=i ^^^IV'pP with Z = Ylit=i e^^- Since there is a uniform 

— Ep 

gap with respect to e between Ei and Ep for p > 2 (see Lemma [4. 2p . the terms e^^ with 
p > 2 are expected to be negligible when compared to the first one {p = 1). Therefore, it 
is natural to expect the solution of ([T]) to be asymptotically close to the solution of the 
following Schrodinger-Poisson system in which only the first energy level is taken into 
account: 






+ V(pi = Siifi 



ZE[0A], 



Si 



— £ 



inf 

<pe/fi(o,i), 11^11^2=1 



W 



/|2 



V^^' 



(3) 



dz'^ 
V{0) = 0, 



0. 



Moreover, when e goes to zero, we will prove that the electrostatic potential, V^, solution of 
(j3]) converges towards a boundary layer potential with profile, Uq, solution of the following 
half line problem: 



d^ij 



- + UiPi = Ei^Pi, ee[o,+oo[. 



de 
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1.2 Main results 

In this paper, a rigourous analysis and comparison of the systems presented above will 
be provided. Namely, ([T]) and ([3]) are posed on a bounded domain. The one dimensional 
Schrodinger-Poisson system on a bounded interval was studied by Nier in [2^. Each of 
these systems can be reformulated as a minimization problem (see section 2 for details). 
However, the limit problem (jl]) is posed on an unbounded domain. Our first result deals 
with the study of (jl]). We also prove that it can be formulated as a minimization problem. 

Theorem 1.1. Let Jo{.) be the energy functional defined on Hq(R'^) (given by (l22l) ) by 

1 r+oo 

^o(^) = 2| \Ur-Er[U], (5) 

where E^ [U] is the fundamental mode of the Schrodinger operator given by (133 p . The 
limit problem (jlj) has a unique solution {Uq, Ei^q.iPi^q), and Uq satisfies the following 
minimization problem: 

MUo)= inf JoiU). (6) 

The comparison of the systems presented above is established by our second main 
theorem. 

Theorem 1.2. Let V^, V^, and Uq be the potentials satisfying problems ([1]), ([3]), and (jl]), 
respectively. Then the following estimates hold: 

l|K-V;iki(o,i)=C?(e-^) (7) 

and 

V;-Vo(-) =0{e-^), (8) 

£ Vfc/ ^1(0,1) 

where c is a general strictly positive constant independent of e. 

The paper is organized as follows. In the next subsection, we present some remarks 
on the scaling giving model ([1]), and we end this section by fixing some notation and 
definitions. In section [21 we recall the spectral properties of the Schrodinger operator on 
a bounded domain and state the optimization problems corresponding to (jl]) and (j3]) (or 
more precisely to the intermediate systems (ITU1) and ( ITT]) ). Section [3] is devoted to the 
analysis of the limit problem (jl]) posed on the half line (proof of Theorem II. ip . We will 
first study the properties of the fundamental mode of the Schrodinger operator (Propo- 
sition 13. ip . The limit problem leads us to the study of a minimization problem posed 
on an unbounded domain. This will be done by means of the concentration-compactness 
principle introduced by Lions in [T7] . Estimates (]?]) and (jH]) are proved in section JH Some 
comments concerning the Fermi-Dirac statistics, the choice of the boundary conditions, 
and the problems of the multidimensional case are given in section Finally, Appendix IX] 
is devoted to the proof of Lemma 13. 1[ 

First, let us make this remark. 

Remark 1.1. To prove (j7])-(j8]), we use the scaled versions of (jl]) and (j3]) when applying 
the changes of variables (j2]) and 

^i(^) = ^V^i f-) , £l = ^E^, V{z) = \u(-), e=-. (9) 
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Remark that it is natural to expect (TTTl) to 6e c/ose, when e goes to zero, to the limit 
problem (j4j) posed on [0,+oo). 



1.3 Remark on the scaling 

Here we show how the system ([1]) can be obtained by a rescahng of the Schrodinger- 
Poisson system written with the physical dimensional variables. Indeed, let (xp(Z),Ap) 
be the eigenfunctions and the eigenenergies of the one dimensional Schrodinger operator 
(the confinement operator) — ^^ + W with homogeneous Dirichlet data: 



h d^Xp 
2m dZ"^ 



+ Wxp = ApXp, 



(12) 



where h is the Planck constant and m denotes the effective mass of the electrons in the 
crystal. The {Xp)p is an orthonormal basis of L^(0, L). The variable Z belongs to [0, L], 
where L is the typical length of the confinement. Denoting by n the electronic density, 
this can be written 



+00 



n{Z) = Y,np\Xp{Z)\'. 



(13) 
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In this formula, |xp(Z)p is the probability of presence at point Z of an electron in the 
pth state. Using Boltzmann statistics, the occupation factor rip is given by 
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where ks is the Boltzmann constant, T denotes the temperature, and Ng is the surface 
density assumed to be given. With this notation we have J^ n{Z) dZ = Ng, which means 
that the total number of electrons in the interval [0,L] (per unit surface in the two 
remaining spatial directions) is given. The electrostatic potential W and the electron 
density n are coupled through the Poisson equation: 

- ^ - ^n (15) 

with boundary conditions 

W{0) = 0, ^{L) = 0. (16) 

In (IT5l) . the constant q is the elementary electric charge and eo, £r are, respectively, the 
permittivity of the vacuum and the relative permittivity of the material. 
Let us rescale the problem flT2|) - (TT6|) by noticing that 

z = je[0,l], WiZ) = ikBT)v(jY Ap = {kBT)£p, Xp{Z) = ^iPp(j 



(17) 
We assume that ^^^2 is of the same order of the thermal energy [kBT). In order to 
simplify the mathematical presentation, we suppose that 

kBT. (18) 
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By inserting (fT7|) into the system ( TT2|) -( !T6|) . we obtain, after straightforward computation, 
the system ([T]) in which e is related to the scaled Debye length: 



3 /-^^^^ \ kBTeoEr .^„. 

' =[t) ^ ^" = V^^i^' ^''^ 

where A^ = ^ is the average volume density of electrons. 

1.4 Notation and definitions 

We summarize in this subsection the different variables and notation used in this paper. 

• For the Schrodinger-Poisson problems posed on [0, 1], z denotes the space variable, 
V denotes the potential variable, and {S,(p) represents any eigenvalue and the cor- 
responding eigenfunction of the Schrodinger operator. For systems posed on [0, -] 
or on M"*", we use ^, U, and {E,iIj) as variables. The same notation with ~, i.e., 
(y,£,(f) or ([/,£', ■0), is used for the variables of Schrodinger-Poisson systems in 
which only the first eigenstate is taken into account. 

• For any real valued function V G L^(0, L), where L > is given (L = 1 or - here), 
we denote by H[V] the Dirichlet-Schrodinger operator 

d^ 
H[V] = -— r + V{x) {x = zotC here) (20) 

dx^ 



defined on the domain D{H[V]) = H'^{0, L) fl i/g (0, L). In addition, the sequence of 
eigenenergies and eigenfunctions of H\y] will be denoted by (£^p[y],'0p[y])pgN. . We 
give in the next section the main properties satisfied by the functions V i— ?■ -Ep[V^] 

and V I—)- il^p\V] for any p G N*. 

• The potentials satisfying ([1]) and (|3]) are denoted by Ve and V^. In addition, 
{£p^£, fp^e) and (^p,£, 'fp,e), withp G N*, represent the corresponding energy couples of 
H[Vs] and H[V^], respectively. In other words, Sp^^ := EplV^], ipp^^ := ^p[K], ^p,e : = 
i?p[T4], and ipp^e := '0p[K]- Similarly, the solutions of flTOj) and fITT]) will be denoted, 
respectively, by {Ue,Ep^e,'^p,e) and {JJs,Ep^e,'ipp,e)- Finally, we fix ([/q, -Ei,o, V'l.o) to 
denote the solution of the limit problem (^. 

Let us now define some spaces which will be used throughout this paper. 
Definition 1.1. (i) For L > 0, we define 

i7i.0(0, L) = {[/ G H\0, L), U{0) = 0} . (21) 

(ii) The space Hq{R~^) is defined as follows: 

ifi(M+) = {[/ G L2^,(M+), U' G L2(M+), U{0) = 0, and U > O} . (22) 

(iii) For any < L < +oo, we shall denote by Sl the set of normalized functions of 
Hq{0,L) with respect to the L'^-norm 

SL=\^eHl{Q,L), j^ V^' = l}. (23) 

Here -ffo(0, L) is the space of H^ -functions vanishing on and L, and when L = +oo 

//i(R+) = {ije H\R+), V'(o) = 0} . 
2 Schrodinger— Poisson system on a bounded domain 

We begin this part by recalling some basic properties satisfied by the eigenvalues and 
the eigenfunctions of the one dimensional Schrodinger operator fl20|) . These properties 
are standard and can be found in [161 1201 1211 129]. The operator H[V] is self-adjoint, 
is bounded from below, and has compact resolvent. There exists a strictly increasing 
sequence (ii^p[y])p of real numbers tending to +cxd and an orthonormal basis of L^(0, L), 
{'iljp[V])p, such that -^plV] G D{H[V]) and 

H[V]^Pp[V] = Ep[V]^Pp[V]. (24) 

For \^ = 0, we have by a simple calculation 

Ep[0] = ^, ^,[0](x) = y|sin(^). (25) 

The eigenvalues -E'pfV^] are simple and satisfy the following characterization (min-max 
principle) [29] : 

Ep[V]= min max ^^^^^''f^'\ (26) 



where Vp {D{H[V])) is the set of the subspaces of D{H[V]) with dimension equal to p, 
and (., .) denotes the scalar product in L^. In view of the min-max formula fl2^ . one can 
verify that for any p &W, Ep[.] is an increasing function, which means that 

Ep[V] < Ep[W] iiV <W a.e. 

Moreover, we have the Lipschitz property, for any real valued functions V, W in L°°(0, L), 

\Ep[V] - Ep[W]\ <\\V- H^|Uoo(o,L). (27) 

Besides, one can prove the following lemma [20] . 

Lemma 2.1. For any p E W , the maps 

Ep[.] : L\0,L) ^ R, ^p[.] : L°°(0,L) -^ L\0,L) 

are Gateaux differentiable, and their derivatives are given, respectively, by 



dEp[V].W= f \ijp[V]\^Wdx and 
Jo 

*Mnw = E E^v] - E,iv] {£ '^*'*"'^) * 



(2S 



foranyV,WeL^{0,L). 

Using the spectral properties of the Schrodinger operator, one can prove the following 
proposition. For details on the proof see 



Proposition 2.1. The systems fllOl) and fllip are well posed. They are equivalent, re- 
spectively, to the following minimization problems: 



and 



JJU,) = inf JJU) (29) 



JeiUe) = inf Je{U), (30) 



where M^ = -. The energy Junctionals J^ and J^ are given by 

UU)-\r\U'?+eHoJf;^e-'^\ (31) 
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p=i 



and 

UU) = ^j^ \U'?-E,[U]. (32) 



Each one of problems f l29|) and fl30l) admits a unique solution. 

Remark 2.1. One can similarly study the systems ([T]) and ^ and prove that each one 
is equivalent to an optimization problem. 



3 Analysis of the limit problem (HD 

The aim of this part is to study the well-posedness of the hmit problem (^ posed on 
the half line. Namely, this part is concerned with the proof of Theorem 11.11 We begin 
with the study of the fundamental mode, -Ef^[.], of the Schrodinger operator. Its main 
properties are listed in Proposition 13.11 

3.1 Properties of the fundamental mode of the Schrodinger op- 
erator on [0, +oo) 

We begin by defining the fundamental mode. 

Definition 3.1. For any real and positive function U G Lj^^iW'^), the fundamental mode 
of the Schrodinger operator is 

E^[U] = inf M^), (33) 

where for any i/j G Soo (defined by fl23|) ) we have 

JuW= / \^f+ / U^'- (34) 

^0 Jo 

One difficulty due to the unboundedness of the interval [0, +oo) is that E^[.] might 

not be an eigenvalue but only the lower bound of the essential spectrum. The following 

proposition gives some properties of E^ [.] and some sufficient conditions on the potential 

for which E'^[.] is an eigenvalue. 

Proposition 3.1. 1. The map U i— )■ E'^\U] is a continuous, concave, and increasing 
function with values in M+ := [0, +oo] satisfying 

E^[U]<\ims\xpU{i). (35) 

2. IfU e L]^M'^), U>0 such that E'^[U] < liminf^^+oo U{0, then E'{°[U] is reached 
by a unique positive function ipi[U], which means that there exists a unique positive 
function i^ilU] G Soo such that E^[U] = Jui^pilU])- In addition, we have 

rITpoo r+oo 

-^[uw = j^ mmi'wd^ (36) 

for any function W in L'^{R'^), the space of bounded functions with compact support 

on R+. 

3. Let U G L;^Q^(R+) be a positive function such that lim^_>+oo U{0 exists, U < lim^_>+oo U{^), 
and E^[U] = lim^_j._|_oo U{^). Then we have 

J poo 

for any W G Lg°(M+). 



4. Let a be an arbitrary positive constant. Then we have 



(37) 



Remark 3.1. There is quite a difference between the third case of this proposition, where 
E'^[U] = lim^_>._|_oo t/((^), and the second case, which includes E^[U] < lim^^+oo t/(^)- 
This result is natural and can be interpreted as follows. The classically allowed region for 
a particle with energy E is the set A= {^ & [0, +oo); U{^) < E}. In the case E > U{^) 
on [0,+oo), the set A extends to +oo so that there is no bound state, while in the case 
E < liiiig_>+oo U{^) the set A is bounded and E is a bounded state energy. 



Lemma 3.1. Let U G H^(^'^) such that E^[U] < liminf+oot/- 

sequences {i/jn)n of problem (l33l) are relatively compact in L^(]R+). 



Then all minimizing 



This lemma is needed for the proof of the second point of Proposition I3.1[ It is proved 
in Appendix |Al The proof is based on the concentration-compactness principle. 

Proof of Proposition 13.11 1. Remark first that for any positive function U, E^[U] 
exists and belongs to R+. It is easy to check, from the definition of -E'i°[.], that it is a 
continuous, concave, and increasing function. To prove inequality (15^ . let i/j G 5*00 be 
fixed and set i/js = \/6tp{5^) for any real positive 5. Then i/js G 5*00, and since E^[U] 
verifies fl33|l . we have 

Er[U]<MiJ5)- (38) 

Moreover, we have Ju{ips) = ^^ /o^°° |^'(0N^ + /o^°^ ^(|)V'^(Ot^^- Then 



limsup Juiips) < 

<5^0 Jo 



+00 



lim sup U 

5^0 



^V'(0^^<limsup/7(0. 

J C^+oo 



Taking the hmsup5_j.o ^^ (I2HD, one obtains inequality (135|) . 

2. Let {^rDn be a minimizing sequence of E'^\U\\ i.e., V'n ^ 'S'oo for any n G N* and 
Juii^n) — ^n-s>+oo E^[U]. The sequence {ilJn)n is bounded in i/Q(M"'"), there exist a function 
i/j G i/Q(R"'") and a subsequence also denoted (■0„) such that (■0„) converges weakly to ip 
in ilfnfM^), and since Ju{-) is weakly lower semicontinuous (it is strictly convex and lower 



semicontinuous) we have Juii^) < lim inf „_i.+oo Juii^n)- 



Then 



(39) 



Besides, the hypothesis E^[U] < liminf+oo U implies that the sequence ('0„)n is relatively 
compact in L^(R"'") (see Lemma 13.11) . Then, up to an extraction of subsequence, ('0„)„ 
converges strongly to ip in L^(]R+). Since ||'0n|li2™+) = 1, for all n, we have ||'0||^2™+) = 1, 
and then if) belongs to S'oo. Therefore, in view of the definition of E'J^f^] fl5^ . E'^\U] < 
JuW and with fl39|) we have E^[U] = Jui,'^)- Let us now show that E'^\U] is a simple 
eigenvalue and the corresponding eigenfunction has a constant sign. Indeed, let ^i and "02 
be two minimizers of Jui.-) on S'oo, i-e., '0i,'02 £ S'oo, such that E'^\U] = Juii^i) = Ju{'^2), 

/~72 ~ 

and let 6 = 



^ + ^. The function 



belongs to S'oo, and we have 



Ju{(p) = ^JuiA) + 7;Ju{^2) 



tpl'4'2 - ^2^1 

2(j) 



E^[U] 
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Since E^[U] < 7^/(0) (0 e ^oo), we get /+°° \M!20!2^^^ = g, which imphes that V^i and 
■02 are proportional, and so i?f^[[/] is simple. In particular, ^0 and l^'l are two minimizers 
of E^[U]; they are then proportional, and since J^ I'j/'p = 1 we conclude that -0 = il^'l- 
We then choose i/JilU] = 1-01, which is positive. This is the unique positive eigenfunction 
corresponding to E^ [U] . To end the proof of the second point of Proposition \'S.1\ let W 
be a compactly supported bounded function {W G L^(M"'')) and remark that for a small 
real t we have E'^^lU + tW] < E'^^lU] + |t|||iy||oo < liminf^^+oo ^(0- In addition, since 
W G LJf , we have lim inf ^_^+oo ^ = liniinf5_i.4.oo(f/(0 + ^^iO)- Then, for any small 
real t, we have E^[U + tW] < lim^^+oo(^(0 + ^^iO)- Therefore, for all bounded and 
compactly supported functions W and for allt G M small, E^ [U + tW] is an eigenvalue. 
Let 0t be the corresponding positive eigenfunction. We have 

p+oo r+oo r+oD 

E^[U + tW]= \'ipt\'^d^+ iU + tW)\iJt\'^d^>E'{^[U]+t W\iJt\'^dC 

Jo Jo Jo 

Similarly, one has 

r+oo 

E^[U] > E^[U + tW]-t W\'ilJi[U]\'^dt 

Jo 

Then, if t is a small nonnegative real (without loss of generality), one can write 







Besides, since ("0*)* is bounded in i/Q(M+), there exists a positive function ipo G Hq such 
that ipt converges weakly to 0O; when t — t- 0^, in Hj^^iM.'^) and strongly in L^^^(]R"'"). By 
passing to the limit t — > 0^ in 

-i/j'^ + {U + tW)iJt = E^[U + tW]i/jt 

we obtain 

-0^^' + U^Jo = E^[U]iJo in V'{0, +oo). 

Since i/jq is positive, we deduce that i/jq = il>i[U]. Finally, to obtain (15^ we just have to 
take the limit t ^ 0+ of (HO]). 

3. Remark first that, since E^[.] is a nondecreasing real function, we have for t > 

E^[U-t\W\]-E^[U] E^[U + tW] - E^[U] E^[U + t\W\] - E^[U] 
t ~ t - t ■ 

Therefore, it is sufficient to prove -^[[/].W = for VF > and W <0 (the general case 
can be deduced by passing to the limit t — )■ 0"*" in the above inequalities). 

(i) Let W G L~(M+) and W >0. Then we have E^[U] < E^[U + tW]. Besides, by 
(135D, we have E^[U + tW] < \im^^+^{U{0 + tW{0) = lim^^+oo f/(0 = E^[U]. Then, 
for allW >0 in L^, E^[U + tW] = E^[U], and the result is proved in this case. 

(ii) Let W G L^(R"'"), let W < 0, and let (t„)„gN be a sequence decreasing towards 0+. 
The sequence {E^[U + tnW])n is increasing and satisfies E^[U + tnW] < lim+oo U = 
E'^[U] for all n G N. Therefore, either it is stationary in the vicinity of +oo and in that 
case -;^[[/].VF = 0, or it satisfies 

E^[U + tnW] < lim U Vn G N. (41) 

+00 
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In the latter case, E^[U + tnW] is an eigenvalue and there exists a sequence {ipn) £ 'S'oo, 
"^n > 0, such that 

Er[U + tnW] = Ju+t„w{i^n)= inf Ju+t„wW- 

Besides, we have E^[U + UW] > E^[U] + t„, J^°° Wij^d^ and 

^1°° [U + tnW] - E^ [U] '■+°° 



< - / WVudi- (42) 

The sequence {4>n)n being bounded in i/Q(R"''), one can find a positive function ip G 
ilfQ(M"'") and a subsequence of ('0„)„ also denoted by ('0„)„ such that f/^n converges weakly 
to ip in Hl^J^'^) and strongly in Lf^^{E.^). In addition ip satisfies, in the sense of distri- 
butions, 

-ip" + Uip = E';°[U]ip = lim([/)?/;. 

+ 00 

This implies that tp" = (U - lim+oo U)tp < with ^ G i/o(M+). We deduce that ^p = 
a.e., and we get the result by passing to the limit in f l42l) . W being compactly supported. 
4. Let us now verify the identity fl37|) . Since the potential (a^/^) tends to +cxd when 
^ goes to +00, E^[ay/^ is reached by a positive function -0 G S'oo: 

Setting ^ = a^C and "PiC) = va^'0(a^C) fo^^ ^^ arbitrary constant /3, we get 

-^^"(0 + "'"-^ v^^(c) = EnaV^mc). 

By choosing /3 such that —2/3 = 1 + |, so that /3 = ^, we obtain 

-^"(C) + yC^(C) = a^Ena^/Cm), 
which implies that E'^[ay^] = a5E^[^/C,]. The proof of Proposition 13 . 1 1 is achieved. 

3.2 Proof of Theorem [TTI 

In what follows, we will show that i^ admits a unique solution verifying (jlj). Indeed, 
the functional Jo(-) is obviously continuous and strictly convex on //q(R"''). To prove 
the existence of a unique Uq G ifo(R+) satisfying (|6]), it remains to verify that Jo(.) is 
coercive on iiro^(M+). For this let U G -f^o(M+) and write U{^) = J^ U'{t)dt. This implies 
that f/(0 ^ II^'IIl^a/^, and since E^[.] is an increasing function one has E^[U] < 

^rlll^'lU^v^]- Applying daZD with a = \\U'\\l2 we get E^[U] < \\U'\\l,E^[^, and 
finally we have 

Mu) > hu'wi, - E^[y^].\\u'\\i — > +00. 

^ \\U\\ui^+oo 

Let us now prove that Uq is a solution of the limit problem (jl]). Namely, we have to check 
that E^[Uo] is an eigenvalue. To this aim, we first write the Euler-Lagrange equation 
for f/o: 

-Uo=^[Uo]>0. (43) 
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Therefore, Uq is a concave function belonging to i^g(M^). It is thus a continuous, in- 
creasing, and positive function on M"*", and hm+oo Uq exists in M. It now remains to check 
that E^[Uo] < lim_|_oo f/o) which will ensure that E^[Uo] := Eiq is an eigenvalue with 
unique positive eigenfunction i/ji^ G S'oo (see point two of Proposition I3.ip . We proceed 
by contradiction and assume that E^ [Uo 
sition [XT] one obtains -^[Uo] = 0. 
positive function in -f^o(ffi' 



lim+oo Uq- Applying the third point of Propo- 
In view of f H5]) and the fact that Uq is a concave 



we deduce that Uq = and min 



f/e/fi 



JoiU) = 0. But 
a simple rescaling argument shows that Jq takes negative values, and so its minimum is 
negative. To prove this claim, we fix a potential U in i/Q(R"'") such that J„ °° |f/'p = 1, 
lim+oo U = +00 and let ipi G ^oo be the eigenfunction corresponding to E'^[U]. For 5 > 0, 
setting U^{^) = e^U{eS,) and i>l{i) = ^/e^jJl{e(), we have 



He 



+ U^r,{0 = e^Er[UWM 



which implies that Ef[U^] = e'^E^[U]. After straightforward computations, we finally 
obtain 



MU' 



dU' 



d^ 



5 r+oo 



d^-Erm = - 



\U'\'d^-e'E'i^[U] 

r-3 



e'E^p] 1 



2Ef°[U] 
which is negative for e small enough. The proof of Theorem 11.11 is complete 



4 Convergence analysis 

The various models presented in the first section of this work are all well posed. In this 
section, we shall estimate the difference between their solutions in terms of e. Namely, 
we have to prove estimates (jTj) and ([8]). The following lemma will be useful. 

Lemma 4.1. Let (f/o, -Ei^o, V^i,o) be the solution of the limit problem (jl]). There exist 
a,b & M^ independent of e such that for all e small we have 



''1,0||l2(A/e,+oo) 



< ae 



-bM^ 



(44) 



with M, 



Proof We have -i/j'Iq + [/o^i,o = £^i,o^i,o such that ^i,o > 0, ^i,o(0 -^c^+oo (^i,o G 
iifg(M"^)), Ei^Q < lim_|_oo Uq, and Uq increases to its limit at +oo. Then one can find two 
nonnegative constants c and 6 independent of e such that for all e small enough we have 



and ^1,0 (Me 
and 



with 



< ce 



^VSM, 



-<o(0 < -#i,o(0 for e e [M„ +oo[ 

. Let S{^) = ce^^^^ and ip = i/ji^ — S. Then we have 5"' = 5S 
-r{0 + SHO<o, ^e[M„+oo[ (45) 

^(Mg) < and ^(+oo) = 0. 
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By the maximum principle, one deduces that -0 < on [M^, +00 [. Thus 

^i,o(0<ce-^« on[M„+oo[, 
which yields estimate (jH]). D 



We begin by proving the second estimate (jH]) of Theorem II .21 For this we will compare 
(see Proposition H?!]) the potentials Uo and Us solutions of Qj and fl30l) . respectively. This 
will be done thanks to an idea consisting of the reformulation of the problems (jlj) and fITT]) 
as minimization problems whose unknown is the first eigenfunction. This is the subject 
of the following remark. 

Remark 4.1. For G SM^(see Definition II. 1|) . where M^ = -, let us set 

Ae{4>)= \4>'m'di+- / 0^(00^(0 min(e,CKrfC (46) 

Jo '^ Jo Jo 

and for G 6*00 

A(0)=/ \4>\i)?di + - / 0^(00^(0 min(^,CKrfC. (47) 

Jq ^ Jo Jo 

The functional A,, satisfies Ai;{\(f)\) = Ae(0) and the convexity property 



A, JtcPl + (1 - t)02 < M,(0i) + (1 - t)A,(02) 



for t G (0, 1), the inequality being strict if |0i| and |02| are not proportional (these prop- 
erties are also satisfied by Aq). The functionals are obvious weakly lower semicontinuous 
on their domain of definition, in such a way that the minimization problems 



and 



A,(0,)= min A,(0) (48) 

Ao(0o)= , min Ao(0) (49) 

4,(iHl{m.+),\\4,\\^2=i 

have unique positive solutions. The problems fHS]) and (1491) are equivalent, respectively, 
to (1111) and (jl]) . Indeed, the functions 0^ and 0o satisfy 

-<Pl + ?7(0,)0. = A^.0, on [0, M,], 

-00 + ^7(0o)0o = /^o0o on R+, 

where /i^ (respectively, fio) is the Lagrange multiplicator associated with the constraint 
\\(/}\\l'2 = 1 and U{(j)s), U{(J)q) denote, respectively, 

U{cPM)= |0.(C)|'min(e,C)c?C, t/(0o)(O=/ |0o(C)l'min(e, C)c?C- 

Jo Jo 

In addition, since 0^ and 0o are positive and the function K{^, C) = min(^, () is the kernel 
corresponding to the Laplacian in dimension one, we have 

{U{(pe),IJ's,<Pe) = {Ue,Ei^e,i'l,s) and ([/(0o) , /iQ, 0o) = (t^O, £'l,0, V^l.o)- 
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Proposition 4.1. The solutions Uq and Ue of IQ and fl30|) . respectively, verify the fol- 
lowing estimate: 

'^{(Je-U,) =0(e-f), (50) 

L2{0,Af.) 



di 



where c is a strictly positive constant independent of e and M^ = - . This yields estimate 



Proof. We start by comparing y4o('0i,o) and As{i/ji,e)- Let Xe G T'(0, +00) be such that 
X^{^) = 1 on [0,M^ - 1], XeiO = on [M^,+oo[, and < Xe < 1- The function 
Xe-^i,o|(o,Af,) belongs to H^{0,M^), and for e small we have ||x£-^i,o||L2(o,Af,) ¥" 0. Let 
/3e = \\Xe-i>i,o\\L'^o,M,)- Then we have j-^Xe-'^Pifi e Sm, and, with (gl]), f3e = I + 0{e~^). 
Then, in view of Remark 14. H the following inequalities can be straightforwardly justified: 



As{i^i,e) < A, ( -x,.^i,o 1 < ^.(V'l.o) + 0{e--) < Ao(^i,o) + 0{e~-), 

where c is a strictly positive constant independent of e. Besides, we have 

^(V'l.o) <Ao{iJi,e) = A,{i>,^,). 

Here and in what follows, we still denote by f/^i^e the extension of f/^i.e by zero on [M^, +00 [ 
when it is taken as a function on M"*". Consequently, we have 



\Ae{^Pl,e) - Ao{l{ji,o)\ = A,(^i,,) - Ao(^i,o) = 0(6"^). 



(51) 



Furthermore, Aq is uniformly convex on 6*00 and i/ji^ realizes its minimum. Then one can 
find a constant cq > independent of e such that 

11^1,0 - i'l,e\\m{R+) < ColAoilpifl) - Ao(^l,e)|. 



In addition, since y4o(V'i,e) = ^e(V'i,e) and with fl^ . one deduces that 



WAfi-AAmi 



0{e- 



(52) 



The potential Uq — Us satisfies 



de 



(Uo-UM) = l^i,o(Or- lV^i,.(Or on [0,M, 



Then, multiplying this equation by Uq — Ue, one obtains after integration by parts 

2 



d , 



< ^(M,){U„(AQ - (MM,)) 



+ sup (|[/o(0-^.(OI) / {\AA-\^iA')di- 

S,e[Q,Me] Jo 



Me 



Moreover, in view of Remark 14.11 we have, for every ^ G [0, M; 



M, 



00 



{Uo-UM)= I (|Vi,or-|V^i,.r)(C)min(e,C)rfC+ / |^i,ormin(e, O^C 

'0 JMs 

"Mg p+00 

^ex I (i^i,or-i^i,.r)(c)rfc+ / i^i,or(c)c?c 

'0 JMe 



<Me 
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and ^{Me 



d^ 



/;r^(Orf^ = /n^i.opdeThen 






{Uo - lie) 



< M, 



L2(o,Me) 



KU 



(IV'l,! 



iV'i,.r)(c)^c+ / i^i,or(c)^c 



< 2M, 



-'1,0 



-^i„ 



e\\L^{0,M^) 



+ 



^'i-,0\\L'2{Me,+oo) ' 



and with fl52|) and (jH]) one obtains fl50|) . Moreover, with the change of variable V^(.) 



^<^s\-j one deduces that 



d 



K 



and then estimate 



\\ holds. 



L2{0,1) 



d 
^ 



(^. - Uo 



L2(o,M^) 



C(e- 



D 



Let us now give the following result, which shows the existence of a uniform gap 
between the first eigenvalue Ei^^ := Ei[Us] and the others Ep[Us]. 

Lemma 4.2. There exists a constant G > 0, independent of e, such that 

Ep[U,]-Ei^,>G Vp>2. (53) 

Proof. Since (£'p[f7e])p>i is an increasing sequence, it is sufficient to show (1^ only for 
p = 2. We argue by contradiction and suppose that |-E'2[f4] — -E-i.el — ?■ as e goes to zero. 
In view of Remark l4.ll we have Ei^^ = Ae{i^i,e) and -Ei^o = ^o('0i,o)- Then, with f lSTj) . 
\Ei^^ — Eifil = 0{e~~). We deduce that E2[Ue] and Ei^^ converge to Ei^q when e — )■ 0. 
Then the eigenf unctions ip2[Ue] and ^/^i^^, corresponding, respectively, to £'2[f^e] and Ei^^, 
prolonged by zero on [i,+oo), are bounded in i/Q(R"'") with respect to e. There exists 
-01 (respectively, 1/J2) G i/o(R"'") such that 'i/'i^e (respectively, "02 [f4]) converges weakly in 
Hl^^{R'^) to ipi (respectively, ?/'2)- By passing to the limit e — )■ 0"*" in D'(0,+cxd), in the 
equations 



■mUe])" + Uei'2[Ue. 






one deduces that 'ipi and 1^2 are two eigenfunctions corresponding to Ei^. In addition, we 
have 

f+00 

|2 



\Juo{i^l,e) -^1,1 



li^l 



+ 







[/o|^l,er-^l,0 



M, 



'0 

< \El,e 





Afe 



l^Ur+ / ^.i^mI -i?i,o + 



(f/o-^.)!^; 



1,£| 



E1.0I+ sup(|t/o-f/.|) 

[0,Af,] 



-^0. 



Then {iPi,e) (and similarly (V'2[f4])) is a minimizing sequence of "-Ei,o = infj/^es^ Juo{ip)" 
Moreover, since £"1^0 < hm+oo Uq (see the proof of Theorem ILII) and applying Lemma 
13.11 (^i,e) and {ip2[Ue]) (up to extraction of subsequences) converge strongly in L^(R"'"). 
Thus, since ■j/'i.e and V'2[f4] are two normalized and orthogonal functions in L^(]R"^) for any 
£ > 0, we deduce that their limits when e ^ 0, ipi and 1^2, which are two eigenfunctions 
of -E'1,0, are also normalized and orthogonal in L^(M"'"). This contradicts the fact that Ei^q 
is a simple eigenvalue. D 
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Proposition 4.2. The potentials Us and Us solutions of flTO|) and (]TT|) verify that 

\\Us-Us\\HH0M.) = O{e-i^), (54) 

where c is a strictly positive constant independent of e. This gives estimate ([7]). 

Proof. Recall first that Us and Us verify, respectively, ( 12^ and fl5Ul) . To prove estimate 
(1M|) . it is sufficient to compare the energies JsiUs) and JsiUs) because we have 

\\Ue - ^.|Ihi(0,M,) < Co\Js{Us) - Js{Ue)l (55) 

where cq is independent of e. A straightforward comparison gives the following inequali- 
ties: 

JsitJs) < UUe) < JsiUs) < JsiUs). (56) 

Besides, we have Us > 0, and Ep[.] is an increasing function; then Ep[Us] > Ep[0] = e^p^Tc"^ . 
Moreover, since Ei^s converges to Eiq., which is then finite, there exists a constant Ci > 
independent of e such that 

Ep[Us]~E,^s>e^p''r^^-Ci. (57) 

Combining (153|) and (1571) , one finds C2 > and C3 > independent of e such that 

Ep[Us] - Ke > C2eVn^ + Cg Vp > 2. 

This implies that 

Bp[C/sl--El[i/e] 



and since 



e ? =(9(e-^), 

p>2 



JsiUs) = JsiUs) + 5^ log 1 + 5^ e P 



we obtain 



p>2 

JsiUs) = JsiUs) + Oie^e-^), 
which leads to (EH) in view of (l55l) and (l56l). D 



5 Comments 

5.1 Fermi— Dirac statistics 

It is more natural to consider Fermi-Dirac statistics in the high density limit (e — )■ 0). 
Here we give some remarks and elements on the limit in this case. The scaled occupation 
factor of the pth state with Fermi-Dirac statistics is given by 

TT-p = fpoiSp — £f), 

where Sp is the Fermi level and fpr, is the Fermi-Dirac distribution 

fFDiu) = logil + e-^), (58) 
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The scaled Boltzmann distribution function, however, is given by /b(u) = e 
Poisson equation in model ([T]) can be written as follows: 



The 






+00 



/ , /-B(FD) {£p- £f)\Vv\'' 
p=l 



under the following constraint on the Fermi energy: 

+00 ^ 

Y,fB(FD){£p-£F) = -,. (59) 

p=l 

In the Boltzmann case, one can explicitly solve (IM?]) with respect to £p, and we have 



e' ES ^-''' ' 



which yields ([T]). The first remark we give in the Fermi-Dirac case is that e^^ cannot be 
expressed explicitly in terms of e"^*". The analysis of the limit can, however, be extended 
to this case but with technical complications that we have avoided in the Boltzmann 
statistics case. When applying the change of variables ([2]) and Ep = \iF the intermediate 
problem flTOj) becomes in the Fermi-Dirac statistics case 



' de 



+ U^p = EpiP, 



prp: 



Ce 



0, 



AeHUo 



fu 

le 



e 



+00 



^p(O) = 0, ^p 



0, 



^plp, 



p^q 



Jpqi 



p=i 



+00 



2 
PI ) 



p=\ 



Ep — ep 



(60) 



.(0,^0, 1(1 



0. 



Since fpo is a regular, positive, and decreasing function on M, the Schrodinger-Poisson 
system in a bounded domain in the Fermi-Dirac case is well posed and can also be 
expressed as an optimization problem; see the work of Nier [20] in the unidimensional 
case and [22] in higher dimensions. More precisely, (!60|) is equivalent to 



where 



7/ I 2 



f/r-r 



p=i 



Je{Ue)= inf UU), 

c/efl'i'0(o,i) 



;,,|Ma^i£M\^|^,_,, - 



(61) 



Ep[U]-ep[U] 
71 



fFD{u)du 



(62) 

Replacing fpoi-) by /b(m) = e ", j£(.) is nothing else but the functional (EI]) modulo a 
constant independent of the variable U. The uniform gap showed in Lemma 14.21 remains 
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correct. Then, as in the proof of Proposition \A.2\ there are two constants ci > and 
C2 > such that 

{Ep - ep) - {El - ep) > CieVvr^ + Ca Vp > 2. 

Since fpD is a decreasing function and log(l + w) ~ -u when m — )■ O"*", then for all p > 2 

/ Ep — ep'\ / Co 2 2 ^1'~^f\ <^2 2 2 ^1~^F 

log(l + e — ^) <log(l + e-^e-"^P" e — ^) < ce'^e-^'P "" e — ^ 



<ce-l*e-'^^p'"'log('l + e" 
where c > is a general constant independent of e. This implies that 

Ep>2 log (l + e" 



log ( 1 + e 




e e^ . 



Thus, a formal analysis shows that, asymptotically when e — t- 0, ( 160|1 is close to a 
Schrodinger-Poisson system with only the first energy level. However, the rigorous anal- 
ysis of the limit, e — >■ 0, of fl6T|) - fl62l) is more technically complicated than the Boltzmann 
case for which the functional J^ has an explicit expression given by f lSTj) . 

5.2 Boundary conditions and higher dimension 

The choice of Neumann boundary condition at 2; = 1 can be justified for modulation dop- 
ing devices (see p]) for which 2; = 1 is in the bulk of the semiconductor and the hypothesis 
of a vanishing electric field is justified. This hypothesis also makes the analysis simple 
because the boundary layer in the limit e — ;■ 0^ is located at 2 = 0. If V^ satisfies Dirichlet 
boundary conditions, then another boundary layer takes place at z = 1. The analysis 
can probably be extended to this case, but the first eigenvalue will have asymptotically a 
multiplicity 2. The multidimensional problem is more complicated, where the location of 
the electrons in the boundary layer may depend on the geometry of the boundary. Such 
problems have been noticed for the Schrodinger equation with a magnetic field by [71 [12] 
and are beyond the scope of our work. 



A Proof of Lemma 13.1 



This appendix is devoted to the proof of Lemma 13.11 We will use the concentration- 
compactness principle. This principle is a general method introduced by Lions [17] to 
solve various minimizing problems posed on unbounded domains. It is shown that all 
minimizing sequences are relatively compact if and only if some strict subadditivity in- 
equalities hold. The proof is based upon a lemma called the concentration-compactness 
lemma. For more details on the principle, we refer the reader to ^7\. Let us begin by 
recalling the concentration-compactness lemma. 

Lemma A.l (concentration-compactness lemma). 1. Let (p„)„>i be a sequence in L^ 
satisfying /?„ > m M and J^pndx = A for a fixed A > 0. Then there exists a sub 
sequence {pnk)k>i satisfying one of the three following possibilities: 
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(i) (Compactness): there exists i/fc G M such that 



Ve > 0, 3R< +00, 

where Br = {x E^\ \x\ < R}. 
(ii) (Vanishing): 



Pn^{x)dx > X- e, 



Vk+BR 



lim sup / pnf.{x)dx = Q. 

k^+00 j^gK Jy+Bn 

(iii) (Dichotomy): There exists a G ]0, A[ such that for all e > 0, there exist ko > 1 
and pI, pI G -^+(IR) such that for k > ko 



\pn^ -(Pfc+Pfc)IUi <^> 



/ p\dx — a 

Jr 



<e, 



p\dx — {\ — a] 



<e, 



where pi has compact support and dist(supp(p^), supp(p^)) -^k +c>o. 

2. If pn = \un\'^ with Un boundcd in H^{W), there exists a subsequence (pn^) such that 
either compactness (i), vanishing (ii), or dichotomy (iii) occurs as follows: there 
exists a G ]0, A[ such that for all e > there exist ko > 1, ul,ul bounded in H^{M.) 
satisfying for k > ko 



( \Wn,-{ul + ul)\\L2<6{e)—^0, 



\ul\'^dx — a < e. 
dist (supp(4),supp(M^)) 



\ul\'^dx — (A — a) 
-)■ +00, 



<e, 



fc— >+oo 
1|2 |v7„.2|2n 



liminf / {|Vw„,r - iVM^r - \^ui\^}dx > 0. 



First, we need to give some notation. For V G i/g(R"'") and e G M"*", we define 



f r+00 

Is = miUv{^), v^Gi/o'(K+), J 



if =e 



(63) 



where Jv{'^) = /q °^ \v'? + /o °^ V^"^, and 






ip =e 



(64) 



where J^{^) = j^^ \^'\^ + V^ j^^ ^'^ and V^ = hminf+00 V. 



Lemma A. 2. Let V G iJo(]R+) such that E'^[V] < hminf+00 V- Then the following strict 
subadditivity inequality holds: 



L<Ia + IZ.^ VO<a<£. 



'a I -^e-a 



(65) 
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Proof. Take v? G H^(M+) such that /^^"^ v'^ 



V^ 



e and let i/j = ^^. Then ^ e //qH^^). 
Jq ip"^ = a, and Jvi'4') = ^Jv{^)- This imphes that £/q, < a/^ for any arbitrary a > 
and e > 0. We also deduce that ela = dh (and similarly el"^ = al^) for any e,a > 0. 
In particular, if a = 1, eli = I^ (or el^ = I^) for any e > 0. Moreover, by definition 
we have h = E^[V] < V^. Then, for all tp e -f/o(M+) such that J^°^ ^^ = 1, we have 
Ji < Jq Iv^'P + V°°. This implies that Ji < I^, and by multiplying by {e — a), which is 
positive ii < a < £, one obtains e/i — ali < (^ — (y)Ii° and inequality f l65|) holds. D 



Proof of Lemma [3.11 Applying the concentration-compactness lemma for (p„)„: p„(x) = 
|'0„(x)p on M"*" and zero elsewhere, there exists a subsequence {pnk)k satisfying one of the 
three cases given by Lemma [A. II If vanishing (ii) occurs, i.e., if 



lim sup / pn^.{x)dx = Vi? > 0, 



which implies that 



lim sup 



y+R 



\iJndx)\^dx = Vi?>0, 



then, for all £ > small enough, one can find a sequence {Rk)k of increasing positive real 

rB 
/O 



such that Jg '"' \^pn^.{x)\'^dx < e for all k. We have 



+ 00 P+OO 

\'^nfdx+ / V^ldx 











oo rRk /■+00 

\^p'J^dx+ / V£jx+ / Vt/j^^dx 

~'0 JRk 

+ 00 



> / l^'J'dx - \\V\Ue + {V^ - e) 

Jo JRi 



i^nJx. 



rik 



'0 JRk 

This implies that there exists S{e), tending to zero when e — )■ 0, such that 

M^nk) > JooiAk) - Ke) > IT - Ke)- 

Now let k go to +00 and e to zero. Then we obtain 

h > /r, 

which contradicts the strict subadditivity inequality ( 1^ . 

Now we assume that {pnk)k verifies the dichotomy case; i.e., there exists a G ]0, 1[ such 
that for all £ > there exist ko > 1, "0^, "0^ bounded in if ^(M"^) satisfying for k > ko 



"rik 



{^l + ^l)\W<5{e) 



e-s>0 



->o, 



\'tl)\{x)\^ dx — a 



<e. 



+ 00 



\i,l{x)\^dx-{l-a] 



<e, 



dist (supp(?/^^),supp(V^^)) 



fe— >+oo 



-)■ +00, 



liminf / {|V^„J2 - IV^^P - |V^^p}rfx > 0. 
'' Jr+ 



7fc 
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One can write (see [T^) 

^nk = -^l + i^k + ^k, where iplipl = ipli^^ = ipli^^ = {) a.e., 

and without loss of generahty, we suppose that supp('0^) C [Rk, +c>o[, where Rk tends to 
+CXD with k. This imphes that 

+00 /"+00 r+oo p+oo 

V^Pl= V\^Pl\'+ V\^IjI\'+ V\cpk\' 

-'0 -'0 -'0 

r+oo r+00 

> / V\^Pl\' + {V- - e) m'-\\V\U6{e) 

Jo Jo 

and 

Jo Jo Jo Jo 

Hence, 

Mi^nJ >7k + M^l) + JMl) - Ke)- (66) 

Besides, let Q^fc = /q °° |V^a;(^)P^^' l^k = Jq°° \4'ki^)\'^^^- For all fixed e > 0, the sequences 
(afc)fc and (/3a:)a: are bounded in R"*". There are subsequences, still denoted by {ak)k and 
Wk)k, which converge in IR+ to a^ and /3s, respectively, where a^ and /?£ belong to IR+ 
such that 

|tte — a| < ^ and |/3e — (1 — «)| < e. (67) 

Inequality flM]) yields 

Jy(^nJ>7fc + /«.+/&-'^(^)- 

Taking the liminffc of the last inequality and letting e tend to zero, we obtain in view of 
fl67|) and the fact that lim inf ^ 7^ > 

which contradicts the strict subadditivity inequality fl65|) . 

Consequently, the sequence (p„j.)fc verifies the compactness case of the concentration- 
compactness lemma which yields straightforwardly that the minimizing sequence (^„)ri is 
relatively compact in L^(]R"'"). 
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